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Abstract 

In order to extract the Cabibbo-Kobayashi-Maskawa matrix element \V u b\ from 
B — > X u lV[ decays, the overwhelming background from B — » X c lVi decays must be 
reduced by appropriate acceptance cuts. We study the non-perturbative effects due 
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Comparisons with recent data from the CLEO, BABAR, and BELLE collaborations 
favor the shape-function approach. 
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1 Introduction 



To test the predictions of the Standard Model for the simultaneous violation of the charge 
conjugation and parity (CP) symmetries in B-meson decays, it is very important to know 
the matrix element \V u b\ of the Cabibbo-Kobayashi-Maskawa quark- mixing matrix [1] 
very accurately. The uncertainties in existing measurements, by the CLEO [2], BABAR 
[31 HI El [6], and BELLE [H [8j [9] collaborations, are dominantly due to uncertainties in 
the theoretical calculation of partial decay rates to be compared with the experimental 
measurements. Experimentally, the inclusive rate Ar u z„(A$) of B — > X u lT>i decays in a 
restricted region A$ of phase space is measured, where the dominant charm background 
is suppressed and theoretical uncertainties are reduced. The theoretical factor i?(A$) 
directly relates the inclusive rate to \V u b\ without extrapolation to the full phase space, 
as 

|T4b| " i?(A$) • (1) 

The uncertainties in the calculation of i?(A$) dominantly originate from the modeling 
of the Fermi motion of the b quark inside the B meson. Most of the recent analyses 
towards the determination of \V u b\ from measurements of Ar u ^(A<3>) [U [5j [8j [9] rely on 
the calculation of i?(A$) by Lange et al. [10]. They use the so-called shape-function 
(SF) scheme, which is an extended version of the original SF approach [ITJ [T2] with many 
effects due to renormalization-group-improved perturbation theory, higher-order power 
corrections from subleading SF terms, etc. But there are many more approaches known 
for describing the non-perturbative B —>■ b transition. We mention the Altarelli-Cabibbo- 
Corbo-Maiani-Martinelli (ACCMM) model [13], one of the oldest models to describe the 
motion of the b quark inside the B meson. In this model, it is assumed that the B meson 
consists of the b quark and a spectator quark, with definite mass m spec and momentum 
Pspec, which is considered quasi-free. The b quark is treated as a virtual particle with 
a mass depending on p spec . Another popular model for describing the non-perturbative 
B ^ b transition is the model of Bareiss, Jin, Palmer, and Paschos based on the parton 
model approach in the light-cone (LC) limit [TH [T5] . All these models, including the 
SF models, contain phenomenological functions of the respective variables describing the 
motion of the b quark inside the B meson with parameters fitted to the 6-quark mass and 
one or two characteristic moments of these functions. Another approach tries to avoid 
these non-perturbative functions by assuming that the fragmentation of the B meson into 
the b quark and the spectator quark can be described as a radiation process off the b 
quark with a proper coupling inserted in the standard soft-gluon resummation formula 
|16j . For a similar approach, referred to as dressed-gluon exponentiation (DGE) in the 
literature, see Ref. [T7] . 

Given the variety of approaches for treating the non-perturbative transition, it is de- 
sirable to make an attempt to compare these approaches with respect to their predictions 
for i?(A$) and other physical observables. In this work, we shall make such a comparison 
between the simple SF approach and the parton model approach in the LC limit, which 
we shall refer to as the LC approach in the following. Such a comparison of the parton 
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model and the ACCMM model has already been done some time ago in Ref. |18j . 

The outline of this work is as follows. In Sect. [2j we give a short introduction to the 
SF and LC approaches. Section [3] contains the results for i?(A$) for three choices of A$ 
underlying recent experimental measurements by BABAR and BELLE. In addition to 
i?(A$), we also present in Sect. [3] distributions in several kinematical variables and com- 
pare them with measured differential decay distributions. Section H] contains a summary 
and the conclusions. 



2 Theoretical ingredients 

2.1 Perturbative differential decay rate 

The differential decay width of B — > X U \V\ has been calculated up to first order in the 
strong-coupling constant a s by De Fazio and Neubert [19] using a fictitious gluon mass 
to regulate soft and collinear gluon contributions. This result has been confirmed using 
dimensional regularization for the soft and collinear singularities in Ref. [20] and by us. 
The quantity of interest is the triple differential decay rate d 3 T / (dx dz dp 2 ) of 

b{ Pb ) -> X u (p) + l(p{) + (2) 

where X u = u or X u = u + g in the case of single-gluon emission and the assigned four- 
momenta are displayed in parentheses. Introducing p = p u + p g and q — pi + p u , we have 
Pb = p + q. The variables x, z, and p 2 are defined as 

_ 2p b ■ pi _ 2p b -p „ 2 _f_ f o\ 

X „ , Z „ , p „ , yo) 

m£ mi mi 

and take the values 

< x < 1, x < z < 1 + x, max(0, z — 1) < p 2 < x(z — x), (4) 

where x — 1 — x. The variable p 2 measures the invariant mass square of the hadronic 
system X u in units of mf, while, in the 6-quark rest frame, x and z correspond to the 
energies of I and X u in units of m b /2, respectively. For fixed values of z and p 2 , x varies 
in the range 

z - ^z 2 - Ap 2 . z + y/z 2 - Ap 1 

< x < . (5) 

2 - - 2 v ' 

Doubly and singly differential decay distributions are obtained by appropriately integrat- 
ing over d 3 T / (dx dz dp 2 ) . The simplest distribution is the spectrum in x, which reads 

P21EI]: 



l^=2x 2 (3-2x) 
I q dx 



-G{x) 



2tt 



where Cf = 4/3, 

G 2 F \V ub \ 2 ml 

1927T 3 ' 



(6) 
(7) 
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with Gf being Fermi's constant, is the total decay rate at leading order (LO) and 



G{x) = ln 2 (l -x)+2 U 2 (x) + -7v 2 



2 9 82 - 153a; + 86a; 
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3 12z(3 - 2z) 



41 - 36x + 42a; 2 - 16a; 3 , . 

+ — — H1 ~ x) ' (8) 

with Li 2 being the Spence function. By integrating over x, one obtains the well-known 
0(a s ) formula for the total decay rate of b — > X u Wf. 



C F a s / 2 _ 25\ 
2?r \ 4 / 



(9) 



Formulas for other doubly differential distributions like d 2 T/(dzdp 2 ) and d 2 T/(dxdz) or 
singly differential distributions like dT /dz and dT /dp 2 may be found in Ref. [T9]. From 
d 2 T/(dzdp 2 ), also the distribution in the hadronic invariant mass Mx can be calculated. 
In the heavy-quark limit, where Pb = (M^/m;,)^, one has 

M\ = p 2 ml + zm b A + A 2 , (10) 

where A = M B — m b . 
2.2 SF approach 

In kinematic regions close to the phase space boundaries, the perturbative spectra are 
infrared sensitive and expected to receive large non-perturbative corrections. Such cor- 
rections are due to the motion of the b quark inside the B meson and are usually referred 
to as Fermi-motion corrections [13]. In the singly differential spectra, such regions are 
1 — x = C(AQ CD /m b ) for the charged-lepton energy spectrum, 1 — z = (9(AQ CD /m fe ) 
for the hadronic energy spectrum, and the low-hadronic-mass region M\ = (9(Aqcd^&), 
where Aqcd ~ 0.5 GeV is the asymptotic scale parameter of QCD. 

One popular method to incorporate Fermi-motion effects is the introduction of a SF 
F(k+), which is supposed to describe light-cone momentum distribution of the b quark 
inside the B meson [TTJ [T2]. The component k + of the 6-quark light-cone momentum 
varies between — m& and A with a distribution centered around k + = and having a 
characteristic width of O (^Aj . The physical 5-meson decay distributions are calculated 
from a convolution of the perturbative 6-quark decay spectra with F(k + ). This is done 
by replacing the 5-quark mass by the momentum-dependent mass m b + k + . Similarly, the 
parameter A is replaced by A — k + (TTJ . Introducing q + = A — k + , the charged-lepton 
energy distribution, for example, is modified to become [19] 



dT .(B XJP.) = 2 ( M - 2E ' dq+ ^J (11) 



dEi Jo Mb — q+ dx 

where dT /dx is the perturbative spectrum given in Eq. (jHJ), x q = 2E\/(Mb — Q+), and the 
charged-lepton energy E\ varies in the range < E\ < Mb/2. The analogous formulas 
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for the distributions in the total hadronic energy and the hadronic mass may be found in 
Ref. [19j and will not be repeated here. Since we wish to calculate the fractional decay 
rate with cuts on E\ and Mx, we need the doubly differential distribution d 2 T/(dEi dMx)- 
This and the triply differential distribution d 2 T j{dE\dMx dq 2 ) are derived analogously 
to Eq. (1111) . After the implementation of the SF, the kinematic variables take values in 
the entire phase space determined by hadron kinematics. For example, the maximum 
lepton energy is ^ max = Mb/2, whereas it is equal to m&/2 for the phase space of the 
perturbative decay rate. 

Several functional forms of F(k+) are available in the literature. They are constrained 
through moments A n = (/c™) of F(k + ), which are related to the forward matrix elements 
of local operators on the light cone [10]. The first three moments are 

4) = 1, A t = 0, ^ 2 = y, (12) 

where /i 2 is the average momentum square of the b quark inside the B meson [22] . In our 
analysis, we adopt the exponential form [23] 

F(k+) = NA~ C (X - A; + ) c e (1+c)fc+/X , (13) 

which obeys A\ = if one neglects terms exponentially small in m&/ 'A. The condition 
Aq = 1 fixes the normalization factor N, and the parameter c is related to the second 
moment as 

X 2 

A* = —■ (14) 
1 + c 

So, the 6-quark mass (or A) and the parameter c (or /j, 2 ) are the two input parameters 
of F(k+). Our choice of A and fi 2 , will be specified in Sect. El when we present our results 
for the cut-dependent partial decay rates i?(A$). 



2.3 LC approach 

Since the B meson is heavy, the momentum transferred in the decay to the final state is, 
in most regions of phase space, much larger than the energy of hadronic binding, which 
is of O(Aqcd)- This suggests that the semileptonic decay of the B meson can be treated 
in a way analogous to deep-inelastic scattering (DIS) in lepton-proton collisions. There, 
LC dynamics dominates DIS and leads to the well-known scaling of the DIS structure 
functions. This is implemented in the parton model, where in LO the structure functions 
are given by the parton distribution functions. These are functions of the scaling variable 
£, which relates the parton four-momentum p q = ^p p to the proton four-momentum p p . 
In an analogous manner, the hadron decay process B — > X u lV\ is modeled by convoluting 
the parton decay process b — > X u lVi with the distribution function /(£) of the momentum 
Pb = £,Pb of the 6-quark inside the B meson according to 

dT(B - XJV t ) = J de/(0 dT(b - X u lv{)\ ph=ipB . (15) 
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This has the consequence that m;, = £Mb is also smeared with the variable £. The 
distribution function /(£) can be expressed in terms of the matrix element of the LC 
bilocal 6-quark operator between 5-meson states as [24] 



/(0 



4vrM| 



d(yp B )e^(B\b(0) 1 -p B (l- l5 )U(0,y)b(y)\B)\ y , =0 , (16) 



where U (0, y) is a gauge link associated with the background gluon field that ensures 
the gauge invariance of /(£). The distribution function /(£) is positive and has non-zero 
values for < £ < 1 only. It fulfills three sum rules [23] . One of them is due to 6-quark 
number conservation and reads 

(17) 



/ de/(o = i- 

Jo 



Reducing the bilocal operator in Eq. (1TB]) to a local one with the help of the operator 
product expansion [23] in heavy-quark effective theory (HEQT), one obtains two more 
sum rules. They determine, up to 0(AQ CD /ml), the mean value \x and the variance a 2 of 
/(£), which characterize the position of the maximum and the width of the distribution 

EE 



Jo 

JO 



M B V 1 + 3^ 



M B 



where 



K h 



1 (^i ^y fe^>. 



2M B 
1 



2M B 
£ 6 = G 6 + K b . 



Ami 

(B\h^¥-h\B), 



2mi 



(19) 



Here, g s = ^4iia s , h is the 6-quark field, G a/ 3 is the field strength tensor of the strong 
force, and D is the covariant derivative involving the gluon field. The matrix elements 
Gb and K b measure the chromomagnetic energy due to the 6-quark spin and the kinetic 
energy of the b quark inside the B meson, respectively. Both are dimensionless HQET 
parameters of C(AQ CD /m^) and are often related to the alternative parameters 



Ai 
A 2 



-2m 2 b K, 
2 
3 



(20) 



The parameter A2 can be extracted from the B*-B mass splitting yielding 
A 2 = (Ml. - Ml) /4 « 0.12 GeV 2 . Values for A x = 



jji 2 ,, introduced earlier, will be 



6 



specified in Sect. |3j when we present our results. If we introduce these two parameters in 
Eq. (fT8"|) . we have 



m b 

M, 



5(A 1 + 3A 2 )~ 



a' 



B 



Ml 



6m 2 



Ai 



'5(Aj 



3ml 



6ml 




(21) 



Of course, the three parameters mb, Ai, and A 2 only constrain the position of the maximum 
and the width of the distribution. For numerical evaluations, one needs the whole function 
/(£), for which we adopt the ansatz [T5], [25] 



(22) 



The parameters a and 6 are determined from the values of and a 2 . The normalization 
factor iV is fixed by Eq. (TPTj) . For 6 = m b /M B and a — >■ 0, Eq. (|22|) becomes a delta 
function, namely /(£) = <5(£ — m b /M B ). In the following, we shall always use Ai and A 2 
as input to determine a and 6 via Eq. (}2~Tj) . 



3 Numerical results 

The large background from Z? — > X c lVi is the main limitation for measuring \V u b\. To 
reject this background, kinematic cuts have to be applied. Depending on these cuts, 
the acceptance for B — ► X u /z7/ decays is reduced. With such acceptance cuts applied, 
the calculation of the B — > X u lV\ decay rate is more complicated and, in particular, 
influenced much more strongly by the modeling of the non-perturbative B —>■ b transition 
than without cuts. 

In recent experimental analyses, four types of cuts have been introduced to separate 
B — > X u lVi decays from the much more abundant B — > X c Zz7; decays. First, various cuts 
on the charged-lepton energy E\ (with or without an additional cut on the invariant mass 
M x of the hadronic system) were used by the CLEO [2J, BABAR [HE], and BELLE [8] 
collaborations. The three other cut scenarios, which were adopted by BABAR [3] and 
BELLE [7J[9] and which we shall consider here, combine cuts on Ei with cuts on Mx, the 
invariant mass square q 2 of the leptonic system [26], and the variable P + = E x — \px\ 
|27j . where Ex and px are the energy and three- momentum of the hadronic system 
X u , respectively. Specifically, they are defined as: (1) E\ > 1 GeV, M x < 1.7 GeV, 
and q 2 > 8 GeV 2 ; (2) E l > 1 GeV and M x < 1.7 GeV; and (3) E t > 1 GeV and 
P + < 0.66 GeV. The corresponding fractional decay rates will be denoted as r\, r 2 , and 
r3, respectively. They all depend on the description of the non-perturbative b — > B 
transition, for which we shall use the SF and LC approaches as discussed in the previous 
section. 
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Both the SF F(k+) and the distribution function /(£) of the LC approach depend 
strongly on the 6-quark mass and much less on the parameters Ai and A2, as we shall see 
below. For these parameters, we choose m b = (4.72±0.08) GeV, Ai = (-0.25±0.10) GeV 2 , 
and A2 = 0.12 GeV 2 . Since the b quark cannot be observed due to confinement, the value 
of nib can only be obtained indirectly from measurements other than that of B — > X u lVi. 
The value of nib depends on the scheme, in which it is defined. For simplicity, we take 
nib to be the pole mass. The scale-invariant 6-quark mass in the modified minimal- 
subtraction (MS) scheme currently quoted by the Particle Data Group [28] as nib = 
m b {nib) = (4.20 ± 0.07) GeV corresponds to m b = (4.78 ± 0.08) GeV at the two-loop 
level. A determination of m b and Ai by fitting B — > X s ^ decay spectra may be found in 
Ref. [29] , with the result that m b = (4.79±Ho) GeV and A x = (-0.2418;?!) GeV 2 . In the 
analysis of their data [9], the BELLE Collaboration used the values m& = 4.60 GeV and 
Ai = —0.20 GeV 2 within the SF scheme. All these values are consistent with our above 
choice for m b and Ai. With these parameters, we calculate the parameters A and c that 
fix the SF F(k+) in Eq. (fl3"j) as well as, via Eq. (j2"Tl) . the parameters a and b that fix 
the distribution function /(£) of the LC approach in Eq. (122]) . In the latter case, we also 
need as input the parameter A2, which we fix as described above. The values of A and c 
in Eq. f[T3"j) and those of a and b in Eq. (|22p are collected in Tables [T] and [21 respectively, 
for m b = 4.64, 4.72, and 4.80 GeV and for Ai = -0.35, -0.25, and -0.15 GeV 2 . 

Table 1: Values of A (in GeV) and c = — 3A/Ai — 1 appearing in Eq. f[T3|) for various 
values of nib (hi GeV) and Ai (in GeV 2 ). 





4.64 4.72 4.80 


\^ A 
Ai^\ 


0.6392 0.5592 0.4792 


-0.35 
-0.25 
-0.15 


2.5021 1.6803 0.9683 
3.9029 2.7525 1.7556 
7.1715 5.2541 3.5927 



Before we can present our results for r\ 1 r^-, and we need to know the change of the 
fully integrated decay rate of B — > X u lVi due to the Fermi motion of the b quark inside 
the B meson. Therefore, we write 

T(B -> XJvi) = r T(b -> XJut), (23) 

where T(b — > X u lv{) is given by Eq. (Q and the deviation of r from unity measures the 
influence of the Fermi motion. The results for r evaluated in the SF and LC approaches 
with the fixed value a s = 0.22 are given in Tables [3] and HI respectively, for the same values 
of nib and Ai as in Tables [T] and [2J We see that, in both approaches, r is approximately 
equal to one. The variation with m b is very small; r mostly depends on Ai. The deviation 
of r from unity is because the factor nib in r [see Eq. (J7J)] is replaced by (nib + k + ) 5 in the 
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Table 2: Values of a and b appearing in Eq. (122]) for various values of mj (in GeV) and 
Ax (in GeV 2 ). 



Ax \ 


4.64 


4.72 


4.80 




-0.35 


0.007950 


0.006940 


0.005895 


a 




0.8941 


0.9094 


0.9245 


b 


-0.25 


0.005911 


0.005215 


0.004493 


a 




0.8861 


0.9014 


0.9166 


b 


-0.15 


0.003604 


0.003212 


0.002804 


a 




0.8780 


0.8934 


0.9087 


b 



SF case and by (C,m b ) 5 in the LC case. It is instructive to approximate these expectation 
values by their lowest non-vanishing moments. In the SF case, we thus obtain for r : 

r « 1 + — 5± 24 

where A 2 is given in Eq. ( Tl4l) . This yields tq = 1.0374 for m b = 4.72 GeV, almost the 
same value as in Table O The derivation comes from the higher moments, which must 
be even smaller. Of course, these results do not imply that the integrated decay rate 
is almost independent of mj. On the contrary, it is proportional to ml and, therefore, 
changes with this factor. Only the influence of the Fermi motion on this decay rate is 
small and feebly depends on m b , as one would expect. Independently varying m b and Ai, 
we have r = 1.0353±g;g^. Table H] exhibits a similar pattern for r in the LC case. For 
our central choice of m b and Ai, it is almost one. It changes very little with m b and more 
with Ai. Over the whole range of m b and Ai, we have r = 1.0044lo!o309- Approximating 
To by the first two non-vanishing moments, we obtain 



25 „ 20 T ^ , 45Ai _ 25A 2 

^~2mJ 



r ^ 1 + y^ + y ir - 1 -6^-2^' (25) 



which yields r ~ 1.0168 for our default values of m b , Ai, and X 2 . Comparison with 
Table H] reveals that, in the LC case, the higher moments are more important than in 
the SF case. Since the error of ro is doubled as compared to the SF case, the error in 
the integrated decay rate is also larger. From Tables [3] and HI we may also conclude that 
parton-hadron duality is realized to good approximation for the total decay rate, r being 
close to unity. 

Next we present our results for the fractional decay rates ri, r 2 , and r%. For the 
SF approach, they are listed in Table [5] for the same choices of m b and Ai as above. 
The central values are r\ = 0.362, r 2 = 0.676, and r% = 0.602. The results for the 
LC approach are given in Table El the central values being r± = 0.360, r 2 = 0.694, 
and r 3 = 0.667. They are similar to the SF case, expect for r 3 , which is larger in the 
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Table 3: Values of r appearing in Eq. (1231) evaluated for various values of m& (in GeV) 
and Ai (in GeV 2 ) in the SF approach. 



Ai ^\ 


4.64 4.72 4.80 


-0.35 
-0.25 
-0.15 


1.0506 1.0484 1.0463 
1.0369 1.0353 1.0338 
1.0225 1.0217 1.0208 



Table 4: Values of r appearing in Eq. (1231) evaluated for various values of m& (in GeV) 
and Ai (in GeV 2 ) in the LC approach. 



Ai \^ 


4.64 4.72 4.80 


-0.35 
-0.25 
-0.15 


1.0350 1.0335 1.0319 
1.0049 1.0044 1.0041 
0.9747 0.9755 0.9759 



LC case. The SF to LC ratios read 1.00, 0.97, and 0.90. Thus, the fractional decay 
rates are remarkably similar in the two approaches and differ only little from the results 
ri = 0.34, r 2 = 0.66, and r 3 = 0.57 obtained in Ref. [10], which were used in Ref. (30] to 
determine \ V u b\ through a global analysis of the available experimental data. As expected, 
the values of r 1; r 2 , and r 3 depend much more strongly on m b than on Ai, both in the 
SF and LC approaches. The variations of with these two parameters are larger in 
the LC approach than in the SF approach. If we express these variations as errors, 
we have n = 0.362^;^, r 2 = 0.676^;^, and r 3 = 0.602±g;?|g in the SF approach 
and n = 0.360j#o29, r 2 = 0.694±g;^, and r 3 = 0.667±[};^ in the LC approach. We 
notice that, in the LC approach, r 3 becomes abnormally small for mj = 4.64 GeV and 
Ai = -0.15 GeV 2 . 

The similarity of the fractional decay rates n, r 2 , and r 3 in the two approaches con- 
sidered here might be related to a fortunate choice of the cut parameters Ef 3333 , M™ ax , 
(g 2 ) min , and p™ ax ; whereas the distributions in Mx, q 2 , and P + for a fixed value of 
Ef 3333 = 1 GeV say, could differ significantly. To elucidate this point, we calculate the par- 
tial decay fractions r±, r 2 , and r 3 as functions of the cut parameters M™ ax , (q 2 ) 333333 , and 
pmax £ or default values of the input parameters, mj = 4.72 GeV (A = 0.5592 GeV) 
and Ai = —0.25 GeV 2 . For this purpose, we define 



n((fl a )-») = f y dq 2 - 



dq 

l r M T* „, dr 



1 

Ei>\ GeV, M x <1.7 GeV 



f 2 (MS? ax ) = - / dM x 
2 V x ' r Jo X dM 



x 



(26) 
(27) 



Ei>l GeV 
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Table 5: Values of ri, r 2 , and r 3 evaluated for various values of m& (in GeV) and Ai (in 
GeV 2 ) in the SF approach. 



Til L 

Ai 


4.64 


4.72 


4.80 






0.3438 


0.3659 


0.3860 


n 


-0.35 


0.6283 


0.6888 


0.7398 


r 2 




0.5411 


0.6223 


0.6908 






0.3386 


0.3617 


0.3824 


n 


-0.25 


0.6082 


0.6763 


0.7330 


r 2 




0.5076 


0.6016 


0.6796 


r-i 




0.3347 


0.3586 


0.3795 


ri 


-0.15 


0.5828 


0.6633 


0.7290 


r 2 




0.4614 


0.5779 


0.6724 


r 3 



1 /- p + ax ,„ dT 



f3 ( p+m 1 = rl + dR 



{2h 

E,>1 GeV 



r Jo dP+ 

which are related to r\, r 2 , and r 3 as 

ri = fi (26 GeV 2 ) - ¥ x (8 GeV 2 
r 2 =f 2 (1.7 GeV), 

r 3 = f 3 (0.66 GeV). (29) 

In Fig. [H fi is plotted as a function of (g 2 ) max for the SF (solid line) and LC (dashed 
line) approaches. We observe that the difference between the two approaches is rather 
small over the whole range of (g 2 ) max , way up to 25 GeV 2 . Later, when we compare with 
experimental measurements, we shall see that the same holds true for the normalized 
distribution (l/T)dT/dq 2 with the above cuts on E\ and Mx- The situation is very similar 
for r 2 , which is shown as a function of M^ ax in Fig. [2j Here, the difference between the 
two approaches is appreciable only for small values of Mj^ ax , for M™ ax <^1.5 GeV. The 
situation is very different for f 3 , which is depicted as a function of P™ ax for the two 
approaches in Fig. [3j We observe that the two distributions coincide at P™ ax « 0.6 GeV, 
where their slopes are very different, however. The result of the LC approach is somewhat 
larger above this value of p™ ax ; way up to P™ ax « 1.2 GeV, while is significantly smaller 
below. As we shall illustrate below, this may be understood by considering the normalized 
P + distribution {l/Y)dY / dP + with the cut E[ > 1 GeV, which is very different for the 
two approaches. It turns out that the choice of E™ m is not responsible for this difference. 

As for measurements of fractional decay rates P(A$), experimental data for the nor- 
malized distributions (l/T)dT /dM x and (l/T)dT/dP + with cuts on E\ have been pub- 
lished and can be compared to the respective distributions evaluated in the SF and LC 
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Table 6: Values of ri, r2, and evaluated for various values of rrib (in GeV) and Ai (in 
GeV 2 ) in the LC approach. 



Til L 

Ai 


4.64 


4.72 


4.80 






0.3476 


0.3674 


0.3856 


n 


-0.35 


0.6201 


0.7305 


0.7878 






0.5668 


0.6997 


0.7657 






0.3395 


0.3601 


0.3788 


n 


-0.25 


0.5578 


0.6942 


0.7743 


r-i 




0.4584 


0.6674 


0.7520 


r-i 




0.3313 


0.3527 


0.3720 


ri 


-0.15 


0.4945 


0.6300 


0.7557 


r 2 




0.1886 


0.6192 


0.7357 


r 3 



approaches (see Ref. [16] for a similar comparison). Specifically, {l/Y)dY / dMx distribu- 
tions with Ei > 1 GeV have been published by BABAR P [6] and BELLE [9]. In Figs. S 
and [5j we compare these measured distributions to our predictions in the SF and LC ap- 
proaches. Both the measured and predicted distributions are normalized to unity in the 
signal region, which is defined by Mx < 2.5 GeV for BABAR [6] and by Mx < 1-7 GeV for 
BELLE |9j. From Figs. H] and El we see that the predictions in the SF approach are in rea- 
sonable agreement with both measurements, whereas the distributions of the LC approach 
are much too narrow and their peaks are much higher than in the measured distributions. 
A similar comparison is performed in Fig. [6]for the normalized distribution (l/T)dT/dP + 
with Ei > 1 GeV measured by BELLE [9J. Both the measured and calculated distribu- 
tions are normalized to unity in the signal region defined by P + < 0.66 GeV. Again, the 
distribution in the SF approach agrees more or less with the experimental data, whereas 
the one in the LC approach is much too narrow. BELLE [9j also presented experimental 
data on the normalized distribution (l/T)dT/dq 2 with E\ > 1 GeV normalized to unity 
in the signal region defined by Mx < 1.7 GeV and q 2 > 8 GeV 2 . These are compared in 
Fig. [7] with the predictions based on the SF and LC approaches. Here, the two theoretical 
distributions are very similar and both agree with the measurement reasonably well. 

Finally, we turn to the charged- lepton energy distribution dT/dEi. In analogy to 
Eqs. ( l26]) - (l28l) . we define the fractional decay rate 

1 r E F** dT 

n(E?n = f l (30) 

In Tables [7] and El we present the values of f 4 (2.3 GeV) evaluated for various values of mj 
and Ai in the SF and LC approaches, respectively. We notice that, for given values of rrib 
and Ai, the results in the two approaches differ appreciably. In particular, r4(2.3 GeV) 
depends much more strongly on rrib for fixed Ai and vice versa in the LC approach as 
compared to the SF approach. This is due to the fact that the Ei distribution falls 
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off much more rapidly towards the threshold at _E' / max = Mb/2 in the LC approach as 
compared to the SF approach (see Figs. l9T fTT]) . Of course, this effect diminishes if Ef 11 ^ 
is taken to be smaller than 2.3 GeV. In this case, also the sensitivity of r± (i^ max ) on m& 
and Ai is reduced. The SF result for f 4 (2.3 GeV) agrees quite well with the value used by 
CLEO [2] to determine \V u b\ from the data points in the range 2.3 GeV < E\ < 2.6 GeV. 

Table 7: Values of f 4(2.3 GeV) defined in Eq. (!30l) evaluated for various values of m b (in 
GeV) and Ai (in GeV 2 ) in the SF approach. 



Ai \^ 


4.64 4.72 4.80 


-0.35 
-0.25 
-0.15 


0.9307 0.9099 0.8874 
0.9419 0.9210 0.8978 
0.9558 0.9347 0.9107 



Table 8: Values of 7^(2. 3 GeV) defined in Eq. ( 1301 evaluated for various values of m& (in 
GeV) and Ai (in GeV 2 ) in the LC approach. 



Ai \^ 


4.64 4.72 4.80 


-0.35 
-0.25 
-0.15 


0.9646 0.9385 0.9106 
0.9780 0.9524 0.9241 
0.9910 0.9669 0.9382 



In Fig. [8], f 4 is displayed as a function of _E™ ax for the SF (solid line) and LC (dashed line) 
approaches. We observe that, as Ef 1 ^ approaches its kinematical upper limit, the LC 
result is saturated appreciably earlier than the SF one. This would lead to an according 
difference in the value of \V u b\ extracted from the data if a large Ef 1 '^ cut were imposed. 
For _E™ ax <^ 2 GeV, the SF and LC results for are very similar. 

In Figs. [91 [TDl and [HI we compare the normalized E\ distributions (l/T)dT/dEi pre- 
dicted by the SF and LC approaches with measurements by CLEO [2J, BABAR [5], and 
BELLE [S], respectively. Both the measured and calculated distributions are normalized 
to unity in the signal region, which is defined by E\ > 2.30 GeV for CLEO [2] and by 
Ei > 2.25 GeV for BABAR [5] and BELLE [8]. In the signal region, where the back- 
ground from b — > c transitions is expected to be minimal, the SF results agree with the 
CLEO, BABAR, and BELLE data quite satisfactorily, while the LC results are clearly 
disfavored. In fact, the E\ distributions of the LC approach drop off much too strongly 
towards the threshold at E\ = Mb/2 and deviate from the data throughout the signal 
region. This disagreement again points to the inadequacy of the LC approach to describe 
the non-perturbative effects in B — > X u lV\ decays, which was already noticed for the Mx 
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and P + distributions in Figs. HH6l Finally, we should note that, in Figs. |9HTT| the theoret- 
ical predictions refer to the rest frame of the B meson, while the experimental data refer 
to that of the T(45) meson. However, since the motion of the B mesons in the T(4S) 
rest frame is non-relativistic, this mismatch is rather insignificant in comparison with the 
experimental errors. 



4 Conclusions 

We studied non-perturbative effects on B — > X u lV\ decays due to the motion of the b 
quark inside the B meson adopting two approaches frequently discussed in the literature, 
namely the shape-function formalism and the parton model in the light-cone limit. While 
these effects are generally small for the total decay rate, they may become substantial once 
kinematic acceptance cuts are applied. In fact, such acceptance cuts are indispensable in 
practice in order to suppress the overwhelming background from B — > X c lVi decays. We 
considered three cut scenarios, involving the invariant mass Mx of the hadronic system 
X u , the variable P + = E x — \px\ related to the energy E x and the three-momentum 
Px of X u , the invariant mass square q 2 of the leptonic system, and the charged- lepton 
energy Ei, that were adopted in recent experimental analyses by the CLEO, BABAR, 
and BELLE collaborations. Comparisons with decay distributions in Mx, P+, and E\ 
measured in these experiments disfavor the light-cone approach. 
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Figure 1: Fractional decay rate r 1 defined in Eq. ( |26|) evaluated as a function of (q 2 ) 
in the SF (solid line) and LC (dashed line) approaches. 
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Figure 2: Fractional decay rate f 2 defined in Eq. (I27p evaluated as a function of M™ ax in 
the SF (solid line) and LC (dashed line) approaches. 
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Figure 4: The decay distribution (l/T)dT/dM x with E x > 1 GeV normalized to unity 
in the signal region (Mx < 2.5 GeV) as predicted in the SF (solid line) and LC (dashed 
line) approaches is compared with BABAR data [6]. 
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Figure 5: The decay distribution (l/T)dT/dMx with E\ > 1 GeV normalized to unity 
in the signal region [Mx < 1.7 GeV) as predicted in the SF (solid line) and LC (dashed 
line) approaches is compared with BELLE data [9]. 
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Figure 6: The decay distribution (l/T)dT/dP + with Ei > 1 GeV normalized to unity in 
the signal region (P + < 0.66 GeV) as predicted in the SF (solid line) and LC (dashed 
line) approaches is compared with BELLE data [9]. 
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Figure 7: The decay distribution (l/T)dT/dq 2 with E\ > 1 GeV normalized to unity in 
the signal region (M x < 1.7 GeV and q 2 > 8 GeV 2 ) as predicted in the SF (solid line) 
and LC (dashed line) approaches is compared with BELLE data [9]. 
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Figure 8: Fractional decay rate f 4 defined in Eq. (1301) evaluated as a function of E™ ax in 
the SF (solid line) and LC (dashed line) approaches. 
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Figure 9: The decay distribution (l/T)dT/dEi normalized to unity in the signal region 
(2.30 GeV < Ei < 2.60 GeV) as predicted in the SF (solid line) and LC (dashed line) 
approaches is compared with CLEO data [2]. 
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Figure 10: The decay distribution (l/T)dT/dEi normalized to unity in the signal region 
(2.25 GeV < E l < 2.60 GeV) as predicted in the SF (solid line) and LC (dashed line) 
approaches is compared with BABAR data [5]. 
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Figure 11: The decay distribution (l/T)dT/dEi normalized to unity in the signal region 
(2.25 GeV < E l < 2.60 GeV) as predicted in the SF (solid line) and LC (dashed line) 
approaches is compared with BELLE data [S]. 
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